Evidence for zero-differential resistance states in electronic bilayers 
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We observe zero-difTerential resistance states at low temperatures and moderate direct currents 
in a bilayer electron system formed by a wide quantum well. Several regions of vanishing resistance 
evolve from the inverted peaks of magneto-intersubband oscillations as the current increases. The 
experiment, supported by a theoretical analysis, suggests that the origin of this phenomenon is based 
on instability of homogeneous current flow under conditions of negative differential resistivity which 
leads to formation of current domains in our sample, similar to the case of single-layer systems. 

PACS numbers: 73.43.Qt, 73.63.Hs, 73.21.-b, 73.40.-c 



Studies of non-linear transport in high-quality two- 
dimensional electron system (2DES) have revealed many 
interesting phenomena which occur in a perpendicular 
magnetic field at large filling factors. In the presence 
of AC excitation by microwaves, there exist microwave- 
induced resistance oscillations (MIROs) [l| which obey 
the periodicity uj/uJc, where w and uJc are the radiation 
frequency and the cyclotron frequency, respectively. The 
minima of these oscillations evolve into zero-resistance 
states (ZRS) for high electron mobility and elevated mi- 
crowave power 0. The MIROs have been found also in 
bilayer and trilayer electron systems [sj, where they in- 
terfere with magneto-intersubband (MIS) oscillations Q 
because of the presence of more than one populated sub- 
band. Recently, it has been demonstrated [5| that ZRS 
exist in bilayer systems despite of additional intersub- 
band scattering. Stimulated by experimental findings, 
theorists have proposed several microscopic mechanisms 
which reasonably explain non-linear transport caused by 
microwave excitation [^-Q. 

A direct current (DC) excitation of high-quality 2DES 
leads to another group of non-linear phenomena caused 
by Landau quantization and, therefore, distinct from 
electron heating effects observed under similar conditions 
in samples with lower mobilities. The Hall-field induced 
resistance oscillations (HIROs) are found in numerous 
experiments [Tol - [l3j and are explained [l^, in terms 
of large-angle elastic scattering between Landau levels 
(LLs) tilted by the Hall field. Further, even a moderate 
DC causes a considerable decrease of the resistance jisj . 
which in high-quality samples may lead to the zero differ- 
ential resistance phenomenon. The zero-differential resis- 
tance states (ZdRS) found in single-layer systems emerge 
either from the inverted maxima of Shubnikov-de Haas 
(SdH) oscillations at relatively hig h magnetic fields [l^ 
or from a minimum of HIROs [131 ■ ^^e second case 
ZdRS appears at low magnetic fields, before the onset of 
SdH oscillations, and extends over a continuous range of 
fields. The two seemingly different regimes, however, are 



explained within the same model assuming formation of 
current domains when the negative resistance conditions 
are reached and homogeneous current picture becomes 
unstable (Tsj . In order to learn more about the origin 
of ZdRS, clear understanding of the domain model is re- 
quired. In this connection, studies of the systems which 
differ from standard 2DES are of crucial interest. 

In this Rapid Communication, we report observation 
of ZdRS in a high-quality bilayer electron system formed 
by a wide quantum well (WQW). Owing to scattering 
of electrons between two populated 2D subbands, ZdRS 
in bilayer electron systems evolve from inverted MIS os- 
cillations for relatively low DC and obey MIS oscilla- 
tion periodicity. In contrast to both regimes of ZdRS in 
single- layer systems [H, [l^l , ZdRS occur in the interme- 
diate regime from overlapping to separated LLs (0.1 T< 
B <0.3 T). A theoretical consideration of magnetoresis- 
tance in DC-driven electronic bilayers shows that nega- 
tive differential resistance (NDR) can be reached in our 
system, thereby supporting the domain model. 

Our experiments are carried out on high-quality 
WQWs with a well width of w=45 nm, high electron 
density Us ~ 9.1 x 10^^ cm~^ and a mobility of /i ~ 
1.9 X 10^ cm^/V s after a brief illumination with a red 
light-emitting diode. Several samples in Hall bar geom- 
etry (length L = 500 and width W = 200 /Ltm) have 
been studied, while we focus here on two samples (A 
and B, with a slightly higher electron density for sam- 
ple B). Measurements have been performed at mK tem- 
peratures in a dilution refridgerator (base temperature 
50 mK) (sample A) and up to 4.2 K in a cryostate with a 
variable temperature insert (sample B). We record longi- 
tudinal resistance using a current of 0.5 fiA at a frequency 
of 13 Hz. Direct current Ijjc was applied simultaneously 
through the same current leads to measure the differen- 
tial resistance r^x = dVxx/ dloc ■ 

In Fig. [1] we show differential resistance Txx as a 
function of the magnetic field for /_dc=0; Idc=^ 
and /£)C'=15 fiA. Dark magnetoresistancc exhibits well- 
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developed MIS oscillations and confirms the existence of 
two populated subbands The inversion of MIS os- 
cillations for i? < 0.17 T is caused by the alternating 
current of 0.5 ^A. The position of the inversion field is 
in a reasonable agreement with our theoretical estimates 
p^ . At this low temperature, we also observe SdH os- 
cillations which are superimposed on MIS oscillations for 
B > 0.15 T. If we apply a constant DC, the MIS oscilla- 
tions are inverted in the whole range of magnetic fields for 
Idc=^ at a temperature of 50 mK. For /£)c=15 /xA, 
we observe ZdRS for B > 0.14 T which are developed 
from the inverted MIS oscillations. Having a closer look 
at magnetic fields 0.1 T< B <0.3 T, we see that the 
maxima between ZdRS are split, which also occurs when 
applying a strong AC or DC excitation separately [l^ 
and can be explained by the influence of Landau quanti- 
zation on inelastic scattering of electrons in two-subband 
systems (see also Fig. SJ. 
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FIG. 1: (Color online) Differential magnetoresistance r^x as a 
function of the magnetic field at Idc=5 fiA and 7dc=15 fiA 
for sample A. The top trace, which is shifted up for clarity, 
shows magnetoresistance {Idc~0) which exhibits MIS oscil- 
lations together with SdH oscillations. Experimental setup is 
shown in the inset. 

As concerns temperature dependence, we have to apply 
a higher DC to maintain ZdRS at higher temperature, see 
Fig m^a) later. The reason for this is rooted in the fact 
that the nonlinear response to a DC bias weakens with 
increasing temperature due to an increase in electron- 
electron scattering d, [13] . 

In order to give a more detailed overview of the evolu- 
tion from MIS oscillations (/_dc=0) to the regime where 
inversion of all MIS peaks takes place until we reach zero- 
differential resistance, we plot in Fig. [l]rj;i^ as a function 
of the magnetic field starting from ljj(j=0 (top trace) 
to J£)(;7=19.2 iiA (bottom trace) at a constant tempera- 
ture of 50 mK for both negative and positive magnetic 
field orientation in steps of AI^c = 0.4 /iA. Starting 
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FIG. 2: (Color online) Evolution of differential magne- 
toresistance r^cx with increasing DC bias from Idc=0 to 
Iijc=W.2 fiA (step AIoc=OA /xA). For the highest current, 
we find four ZdRS for B > 0.12 T. 



from /£)C'=16 ^A the widths of ZdRS remain constant. 
A further increase in Ijjc leads to heating of the electron 
gas in the sample, therefore our studies at r=50 mK are 
limited to a maximal current of I^q^=20 /iA. 

Since ZdRS in electronic bilayers occur neither at 
nearly discrete values of B [l^ nor in a wide region of 
magnetic fields fl7^ , and are separated by maxima in r^^ , 
it is important to illustrate this phenomenon by sweep- 
ing luc at a constant magnetic field, see Fig. [3] We 
fix the magnetic field at several positions (split maxima 
and ZdRS in rxx)i mark them by arrows with numbers 
(1) to (5) in Fig. EJa), and plot differential magnetore- 
sistance r^^ilDc) in Fig. [3i;b). For B = 0.254 T and 
B = 0.189 T, we are situated in the center of ZdRS. Both 
ZdRS are developed at Idc > 10 ^J.A. For B=0.254 T we 
find small kinks in r^x at ±5.8 /iA before r^x is stabilized 
near zero. Beyond ZdRS, we show similar plots for the 
two maxima in Vxx at B = 0.204 T and B = 0.224 T 
which occur between the ZdRS, and in the transition 
from these maxima to ZdRS at i? = 0.2 T. The longitu- 
dinal voltage Vxx as a function of Idc is plotted in Fig. 
M^c) for the ZdRS around B:=0.25 T. This trace is ob- 
tained by integrating the data at B = 0.254 T and shows 
a typical voltage-current characteristic with saturation of 
Vxx corresponding to the onset of ZdRS. 

We now analyze and discuss our experimental observa- 
tion. First, we point out that the occurrence of ZdRS in 
electronic bilayers is distinct from the case of single-layer 
systems where ZdRS emerge in a wide range of magnetic 
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FIG. 3: (Color online) (a) Differential resistance r^x for 
/_DC=20 jj,A. (b) Dependence of r^x on the DC at several 
chosen magnetic fields as labeled at a temperature of 50 mK. 
(c) Corresponding voltage Vxx at B=0.254 T. 



fields well below the onset of SdH oscillations. In a two- 
subband system several ZdRS emerge from the peaks of 
inverted MIS oscillations and follow the 1/B periodic- 
ity of these oscillation determined by the subband sep- 
aration A. In other words, each ZdRS appears around 
the magnetic fields given by the set of cyclotron ener- 
gies hujc = A/fc, where k is an integer. Next, the other 
important properties of ZdRS, such as the dependence 
of r,j:x on the DC and the voltage-current characteristics, 
are very similar to those reported for single-layer systems 

Mm 

From the theoretical point of view, it is quite pos- 
sible that the differential dissipative resistance r^^ = 
dVxx/dlDC (or even absolute resistance Rxx) in high- 
mobility samples approaches zero with increasing cur- 
rent, since an application of the current leads to a con- 
siderable decrease in the resistivity. Such a decrease, 
observed both in single-layer [l5| and double-layer [l^ 
systems, is reasonably explained by the theory of Ref. 
(9| based on a modification of the electron distribution 
function due to current-induced diffusion of electrons in 
the energy space in the presence of Landau quantiza- 
tion. This non-equilibrium distribution function is sta- 
bilized by inelastic electron-electron collisions described 
through the inelastic scattering time (inelastic mech- 



anism). A more rigorous theory, which includes another 
(displacement) mechanism of nonlinear response and re- 
mains valid at higher currents, has been developed in 
Ref. ■ To investigate the possibility of NDR for our 
samples and to explain other characteristic features of the 
observed magnetoresistance, we have generalized the the- 
ory of Ref. |14| beyond the regime of overlapping Landau 
levels and for the case of multisubband occupation. Be- 
low we present the results valid both for single-subband 
systems and for the systems with two closely spaced sub- 
bands (our samples). In the regime of classically strong 
magnetic fields, the magnetoresistance of the system of 
electrons interacting with the static potential of impu- 
rities is found from the following relation between the 
density of the applied current j = Idc/W and longitu- 
dinal electric field E\\ = Vxx/L: 
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where po = rn/e^nsTtr is the classical Drude resistiv- 
ity (e and m are the electron charge and the effec- 
tive mass). The contribution of SdH oscillations is ne- 
glected. The sums are taken over the cyclotron har- 
monics, k are integers, au are the coefficients of expan- 
sion of the density of electron states in such harmonics, 
D{e) = Cfc exp(ifc|^), where D[e) is the density of 
states in units of mjiif? . For two-subband systems, one 
should use D{e) = {Die + D2e)/2, where Di^ is the den- 
sity of states in the subband i. In both cases ak are 
taken real, which is attainable by a proper choice of the 
reference point for energy. The coefficients Aki , which 
describe harmonics of non-equilibrium oscillating correc- 
tion to the distribution function, are found from the lin- 
ear system of equations ^j., Ckk'Aw ~ akl'k with 

Cfcfc' = —[a%_i.,+af.{ak+k'—'iak-k')]+ak-k'{lk-k'^lk)- 

(2) 

The first part of the matrix Ckk' is determined by tak- 
ing into account the Landau quantization in the lin- 
earized electron-electron collision integral (see details in 
Ref. [13 )• The coefficients ^k describe the effect of 
the current and are defined by averaging over the scat- 
tering angle 0: 7fe = 7(Cfc) = Jo[2Cfc sin (g/2)]Tt jT(g), 
where Jo (a;) is the Bessel function, C, = \J Air'^p /e^ngool 
(for single-subband systems one should multiply C by 
\/2), and 1/t{0) = {m/rt')w[2kF sixi{e / 2)] is the angular- 
dependent scattering rate {kp is the Fermi wavenum- 
ber, w{q) is the correlator of random scattering poten- 
tial). The transport time r^^ is given in the usual way, 
1/Ttr = (1 — cos 9) /t{9). Finally, 7^ and 7^' denote, re- 
spectively, the first and the second derivatives of the func- 
tion 7(Cfc) over its argument. 

To carry out the calculations, we numerically com- 
puted the density of states in the self-consistent Born 
approximation, using the quantum lifetime of electrons 
To = 6.6 ps determined from linear low-temperature 
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FIG. 4: (Color online) (a) Differential magnetoresistance r^x 
in experiment and (b) theory for three chosen temperatures 
(sample B). 



magnetoresistance measurements. Next, we applied the 
theoretical @ estimate Tin — fisp/T^, where Ef is the 
Fermi energy and Te is the electron temperature. Heat- 
ing of electron gas by the current is taken into account 
by considering energy loss due to electron-phonon in- 
teraction [2^ (for Idc ~ 50 iiA the heating is ap- 
preciable only at low temperatures, at T = 1.4 K we 
find Te ~ 1.6 K). To describe clastic collisions, we used 
the model of long-range scattering potential, 'w{q) cx 
exp{—lcq) (where Ic ^ l/kp is the correlation length), 
which produces Ttr = to{1 + x) / X with x = ^/{kplcY 



and 7(C) ^ (Tt./To)/yrT^ M- 

In Fig. |4]the calculated differential magnetoresistance 
plots arc compared with experimental data recorded on 
sample B with a subband separation energy A = 1.4 
mcV. Notice that this separation is slightly larger than 
in sample A, so the positions of inverted MIS peaks (and 
of the corresponding ZdRS) are slightly shifted towards 
higher magnetic fields with respect to sample A. Apply- 
ing Idc = 50 /uA, we observe ZdRS around B = 0.27 T 
at low temperatures (see the plot for T = 1.4 K) and 
its disappearance at T > 2.5 K. The calculation shows a 
similar behavior with a NDR region around B = 0.27 T. 
The other important feature of magnetoresistance, the 
peak splitting leading to minima at 0.23 T and 0.32 T, 
is also reproduced in the calculations, see Ref. [l^ for 
a detailed description of this phenomenon. Therefore, 
a direct calculation of magnetoresistance demonstrates 
that NDR is attainable for our samples. In the NDR re- 
gions, according to Rcf. [l^, the homogeneous flow of 
the current becomes unstable and the system exhibits a 
transition to current domains. The simplest two-domain 
structure discussed in Refs. [H, [l^ can as well describe 
ZdRS in our experiment. 

To summarize, we have found evidence for zero- 
differential resistance in a bilayer (two-subband) elec- 
tron system where ZdRS develop from inverted magento- 
intersubband oscillations under a relatively small direct 
current. This experimental result, together with a the- 
oretical consideration, suggests that the ZdRS occurs as 
a result of current instability under negative differential 
resistance conditions. The domain model proposed for 
explanation of ZdRS in single-subband 2DEG [H, [T^ 
very likely applies to multisubband systems. 

Note added in proof. We have recently become aware 
of related experimental work on a similar system, by Ref. 
21. 
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